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Abstract 

After a brief review of the kinematics of deep inelastic scattering 
(DIS) within the so-called £ method, we derive the necessary formulae 
for the treatment of QED radiative corrections to DIS originating 
from hard photon radiation. The results are applied to a calculation 
of the corrections to DIS with a tagged photon with next-to-leading 
logarithmic accuracy under HERA conditions. It turns out that the 
next-to-leading logarithmic corrections are quite important for the £ 
method. We also discuss the dependence of the corrections on the 
longitudinal structure function of the proton, Fl, in the region of low 
Q 2 and moderate x. 
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1 Introduction 



The determination of the structure functions of the proton, F 2 (x, Q 2 ) and 
Fi(x, Q 2 ), over a broad range of the kinematic variables belongs to the most 
important tasks of the HI and ZEUS experiments at the HERA ep collider. 
Especially the extension of these measurements to the range of small Bjorken 
x < 10~ 4 and Q 2 of a few GeV 2 is of particular interest, as it provides a testing 
ground for our attempts to understand the details of the dynamics of quarks 
and gluons inside the nucleon. 

Whilst the structure function F 2 can be extracted quite easily from the ex- 
perimental data, it is more difficult to determine the longitudinal structure 
function Ft- A direct method that relies only on measured data requires 
running the collider at different center-of-mass energies. However, besides 
impairing the high-energy program of the machine, running at reduced beam 
energies also increases some systematic errors, (e.g., luminosity uncertain- 
ties), in the experimental analysis. 

These problems are circumvented by employing a method suggested by 
Krasny et al. M that utilizes radiative events. This method takes advantage 
of a photon detector (PD) in the very forward direction, as seen from the 
incoming lepton (electron or positron) beam. Such a device is part of the 
luminosity monitoring system of both the HI and ZEUS experiments. 

The idea of this method is that emission of photons in a direction close to 
the incoming lepton corresponds to a reduction of the effective beam energy. 
This effective beam energy for each radiative event is determined from the 
energy of the hard photon observed (tagged) in the PD. Early analyses that 
make use of these radiative events for a determination of F 2 were already 
published in |2|,[|. No QED radiative corrections were taken into account in 
these analyses. The feasibility of a determination of Fl was studied in [|]]. 

Recently, the HI collaboration presented preliminary results of a refined 
analysis with newer data 0. In this analysis, the authors chose different 
methods of determination of the kinematic variables^ (the e-method, where 
the kinematic variables are obtained from a measurement of the scattered 
lepton, and the £ method) in different (x, Q 2 ) bins in order to reduce the 
experimental systematic error. However, since the calculations of the QED 
radiative corrections to DIS with a tagged photon |7]-|9| did not cover the £ 
method, the corrections were only applied to part of the data in ||. It is the 
purpose of the present work to extend these analytical calculations to the £ 
method. 

WFor a discussion of the most common methods to determine kinematic variables and 
further references, see e.g. Q. 
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The E method, as proposed by Bassler and Bernardi |nj , tries to combine 
the momenta of the outgoing lepton and hadrons judiciously in order to 
reduce experimental systematic uncertainties on the determination of the 
kinematic variables especially in the kinematic region of low Q 2 where other 
methods are limited by e.g., detector resolution or energy calibration. 

With the help of the quantity!^] 



Pz,h) 



where the sum runs over the detected hadrons, and Eh and p z ^ are the energy 
and ^-component of the respective particle, the kinematic variables x%, Ut, 
and Q| are defined via 

_ S fe _ E' 2 sin 2 9 _ Ql 

yE "£, + ^(i + cos£)' gs "T^' X *~^S- [2} 

Here S = 4E e E p , where E e and E p are the beam energies of the lepton and 
proton beam, respectively, and E' e and 9 are the energy and scattering angle 
of the outgoing lepton, measured with respect to the direction of the initial 
lepton. 

One of the known advantages of the E method is its insensitivity of the 
determination of and Q| to undetected emission of a hard photon collinear 
to the incoming lepton (initial state radiation, ISR). 

The X method has already been used in several analyses of the HI col- 
laboration. However, the author is not aware of any publications on an 
analytical (i.e., non- Monte Carlo) treatment of QED radiative corrections 
to deep inelastic scattering (DIS) using the X kinematic variables^ beyond 
the collinear (leading log) approximation [12, I3|. The presumable reason 



is that the E method has only been introduced long after the start of data 
taking at HERA. Therefore, this paper starts with a brief introduction to the 
kinematics of radiative DIS in the E method. Section ^ extends the consid- 
erations to the case of DIS with an exclusive tagged photon, but specialized 
to the conditions at HERA, and provides the relevant formulae to calculate 



( 2 >Note that in this paper we take the positive z-axis along the initial lepton direction, 
unlike flo|| who chose the direction of the incoming proton beam. 

( 3 )As the tagged photon cross section represents a radiative correction to the DIS cross 
section Q , one can in principle calculate the QED corrections to the former for any choice 
of kinematic variables with the help of a Monte Carlo event generator for DIS that properly 
implements the necessary higher order QED corrections to DIS. However, no generator 
exists for the calculation of QED corrections beyond leading logarithms, and the leading 
log generator KRONOS pd| ] uses approximations for photon emission in the very forward 
direction that make it useless for the present task. 
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the radiative corrections to the tagged photon cross section, based on the 
results of §|. Some results for HERA experimental conditions are presented 
in section |], and section ^ contains our conclusion. Finally, the appendices 
collect several technical details. 



2 Kinematics in the E Method 

This section is devoted to a basic review of the £ method. Here we shall 
prepare an appropriate framework for the treatment of radiative corrections 
to radiative deep inelastic scattering, 

e(p) + P(P) - e{p') + X(P') + 7 (k) , (3) 

i.e., DIS with single hard photon emission. The extension to the process with 
an additional tagged photon in the forward direction is straightforward and 
will be performed in the next section. 

Let us begin by stating our conventions for the kinematics in the HERA 
lab frame that are used throughout this paper. We shall take the orientation 
of the coordinate frame such that the positive z-axis points in the direction 
of the incoming lepton beam, and the momentum of the scattered lepton lies 
in the x-z-plane: 

P = (E p ,0,0,- Pp ), 

p = (E e ,0,0,p e ), 

p' = (E' e ,p' e sin 0,0, pe cos 0) , 

k = E 7 (l, sin i? cos <p, sin $ sin (p, cost?) . (4) 

As the beam energies E p , E e , as well as the energy of the scattered lepton, 
E' e , are always large compared to the proton mass, M, and the electron mass, 
m, we shall take p p = E p , p e = E e , and p' e — E' e , wherever possible. 

Since we assume E p 3> M, we may replace the definition of the variable 
Z h in ©by 

P-(P'-P) 
■= = • (5) 

This allows us to similarly reexpress the definitions (0) of the kinematic 
variables xs, yn, and Q| through scalar products of four- momenta via 

P-(P'-P) 

2/2 = p. { p>- P+pr ^ 

2 4(p.p')(P.p') 

Qs(i-ys) = g , (7) 



3 



^ = % (8) 

where 

S = 2P-p. 

One important thing to note here is the nonlinear dependence of the kine- 
matic variables y% and Q| on the energy and scattering angle of the outgoing 
lepton, while it is linear in the electron-only method (e- method). 

As we are dealing with the kinematics of a process with real photon 
emission, it is convenient to define an invariant quantity k, that represents 
the energy of the outgoing photon in units of the energy of the incoming 
lepton, as measured in the rest frame of the incoming proton, 

P-k 

k := — . 9 
P ■ p 

Energy and momentum conservation of the process ([!]) obviously requires 
< k < 1. In the special case of emission of the photon collinear to the 
incoming lepton, k also represents the energy fraction of the initial lepton 
taken by the photon in the HERA lab frame. 

Using momentum conservation (|3]) and relations (|6|) and (|9|), it is easy to 
see that 

P-(P'-P) = {l-K)yxP-p = (l- K )yx.±S. (10) 

To proceed, let us express the remaining scalar products of external mo- 
menta in terms of invariant kinematic variables and other measured quanti- 
ties. We find 

P . p > = hzMp.( P '-p) = (i_ K )(i_ ys) .I 5) 
, vl S 2 x s y s l 

P - P = AP ■ (P> - P) = —-2 s - 

Furthermore, from the photon's energy E 1 =: x^E e and angles •& and with 
respect to the incoming and final lepton, respectively, we obtain 

p-k = x^E 2 e (1 — cos-#) , 

p'-k = x y E e E' e {l - cos^') , (12) 
where the angle is calculated from (||), 

cos = cos 9 cos i? + sin 9 sin d cos tp . (13) 
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Next, we rewrite the left-hand sides of the relations ( |TT| ) in the HERA 
lab frame as 

P-p' = E p E' e (l + cos6), 

p-p' = E e E' e (l- cos 9), (14) 

to derive explicit expressions for the energy and scattering angle of the out- 
going lepton in the HERA frame, 

E' e = (l- K )(l-yv)E e + ^E p , 

1 — K 

cos 9 = j 1 ~ "ff 1 ~ " X ^ E J , (15) 

that may further be used to eliminate E' e in the second equation of (|T^), or 
to trade the azimuthal angle ip in favor of the angle 

It is now straightforward to obtain the expressions for the momentum 
transfer to the hadronic system 

Ql = -{P-P'f = -{p-p'-kf = ^- + 2k.{p-p'), (16) 

I — K 

the hadronic scaling variable 

Xh = 91 = 91 (i 7 ) 

h ~ 2P-(p-p'-k) (l-K)y E S" 1 1 

and the invariant mass of the hadronic system, 
W 2 = {P'f = (P + p-p'-k) 2 (18) 
= M 2 + y s 



K 



l-« 



S-2k.(p-p') = M 2 + l —^Q 2 h . 

x h 



In the last equation we explicitly retained the proton mass. 

The kinematic limit for the phase space of the radiated photon may be 
derived by requiring that the invariant mass of the hadronic system be larger 
than the threshold for pion production, 

W 2 > M 2 , where M = M + . (19) 

The actual upper limit on the photon energy as a function of the emission 
angle, _E™ ax = E™ ax (9; ip) = E™ ax (9; is obtained in general by solving 
a quadratic or cubic equation; for details see appendices |A] and [FJ. 
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Before we conclude this section, we shall give the relations between the 
kinematic variables in the S method and in the e-only method. From the 
scalar products flTT|), we find: 

1-2/e = (1-k)(1-3/e) =>- 2/e = 2/s + «(1 - 2/s) , 

XeVe ~ 1-k ^ Xe - + { ] 

The Jacobian between these two sets is 

j = det = — ^ — > = PD 

\d(x s ,y s )J yz + K{l-yx) y e 

The reader may easily verify that the above formulae are consistent with 
the collinear limits discussed in Jl2| . 



3 DIS with a Tagged Photon 

After having discussed the kinematics in the E method, let us now turn to our 
primary aim, the description of radiative corrections to neutral current deep 
inelastic scattering with an exclusive tagged photon for the HERA collider. 



3.1 Kinematics and Lowest Order Cross Section 

As already explained in the introduction and described in more detail in [I]] 
(and references cited therein), the experimental detection of photons emitted 
in the very forward direction is actually possible at the HERA collider due 
to the presence of photon detectors (PD) that are part of the luminosity 
monitoring systems of the HI and ZEUS experiments. These PD's cover 
an angular range very close to the direction of the incoming lepton beam, 
(tf 7 = Z(p, k) < $o « 5 • 1(T 4 rad). 

Thus, the process under consideration corresponds to the reaction 

e(p) + p(P) - e{p') + 7 (A0 + X(P') + [ 7 (fe)] , (22) 

where j(k) denotes the collinearly emitted, exclusively measured (tagged) 
photon, and [7(^2)] indicates an additional (i.e., second) photon in the case 
of the radiative correction to the lowest order process. 

A set of kinematic variables that is adapted to take into account tagged 
collinear radiation is given by the shifted Bjorken variables [|1], [7], |9| -0 Ex- 
pressed in terms of the measured quantities, S/j (see eq. ||), the energy E' e 

( 4 )Here we shall use the notation x etc. of |J to avoid confusion between shifted Bjorken 
variables and the usual ones. Although x etc. are determined in the S method, we drop 
the index S in order to not overload the notation. 
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and angle 6 of the scattered lepton, and the energy of the tagged photon, 
they read 



y 



E fc + + cos i 



E' 2 sin 2 9 



x 



yzS 



(23) 



Here we denoted by z the energy fraction of the lepton after initial state 
radiation of a collinear photon, 



2P-(p- k) E e - E n 



E„ 



(24) 



where E n represents the energy deposited in the forward PD. 

The above definition (p3| ) corresponds to eqs. ©-(§), but with an effec- 
tively reduced center-of-mass energy, zS. It is obvious that in the case of 
the £ method only x is affected by collinear initial state radiation, as was 
already found in JlO| . 

In analogy to the case of ordinary DIS, we may rewrite the kinematic 
variables in a Lorentz invariant fashion: 



P-(P'- P) 
P-(P'-P+ p') 



Q 



2 4(zp-p')(P-p') 



(l-y)zS 



x 



yzS 



(25) 



The Born cross section, integrated over the solid angle of the photon 
detector (0 < $ 7 < $ , $ <C 6) takes a factorized form (see also |]il,|7jP,|14]j): 



1 d 3 a Bo 



y dx dy dz 2n 



— P{z,L )Y>. 



(26) 



where 

t = t(x,y,Q 2 ) 

2vrft 2 (-g 2 ) 



Q 2 xy 2 



M 2 



2(l-y)-2x 2 y 2 -^-+ (l + Ax 

Q 



Q2 J 1 + R 



with 



P(z,L 



Q 2 

R = R(x,Q 2 ) 



1 + z z 
1 - z 

xyzS . 



Ln- 



2z 

ft(-g 2 



L = In 

a: 



m 2 



1 + Ax 



2 M 2 \ F 2 {x,Q 2 



Q 2 J 2xF 1 (x,Q 2 ) 
7 



i-n(-g 2 ) 
i . 



(27) 



The quantities F2 and F\ denote the proton structure functions. Note that 
we explicitly include the correction from vacuum polarization n(— Q 2 ) in the 
virtual photon propagator, and that we neglect the contributions from Z- 
boson exchange and 7-Z interference, because we are interested mostly in 
the kinematic region of small momentum transfer Q 2 . 

The cross section fl26|) describes the process (p2|) to lowest order in pertur- 
bation theory. The radiative corrections to this cross section are composed 
of contributions by corrections due to virtual photon exchange, soft photon 
emission, and emission of a second hard photon, with one of the hard photons 
being tagged in the PD. Because of its coarse granularity, we shall assume 
that the PD cannot measure photons individually but only their total energy 
when two hard photons simultaneously hit the PD in different locations. 



3.2 Virtual and Soft Corrections 

The virtual and soft corrections to the lowest order cross section are simply 
obtained from the calculation for a measurement using the e-method by 
substitution. For the sake of completeness, we quote the result from Q: 

with 

A 2 7T 2 9 /l + c\ 

p = 2(L Q -l)ln— + 3L Q + 31n£-ln 2 y---- + 2Li 2 (^-J, 

nn l+Z 2 ,*, „N T , 2-(l-,2) 2 

T = — -(A In Z + B)-— -L Q In z- ^ _ ; L + O (const ) , 
A = - L 2 + 2L L Q - 2L \n(l - z) , B = [in 2 z — 2Li 2 (l — z)] L, 



<1 





t 

O 2 f du 

ln-^r, Li 2 m = - / — ln(l-u). (29) 
" J u 



zm 2 



Here A denotes the infrared cutoff for the emission of a soft photon in addition 
to the hard one, A = E^^/Ee. Furthermore, 

v _ E' e E p _ z(l - y)E e - xyE p 

Y = — — z(l — v) +xy — p and c = cos6 l = — — — - ~r 30 

E e 1 y) y E e z(l-y)E e + xyE p 1 ; 

follow from the formulae of the previous section with the replacement p zp, 
as the energy loss due to the tagged collinear photon is known and already 
taken into account in the determination of the kinematic variables. 
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It should be noted that in (|28|)-([29"|) and also further below we retain only 
terms with double or single large logarithms of the small electron mass m, 
i.e., terms of order a 2 L 2 and a 2 L, with L being one of L or Lq. As the 
lowest order cross section (^) is of order aL relative to the DIS cross section, 
we shall denote the terms of order a 2 L 2 as leading (LL) and those or order 
a 2 L as next-to-leading logarithmic (NLL) ones. 



3.3 Double Hard Bremsstrahlung 



Besides the soft and virtual corrections to the lowest order process, we have to 
consider also the corrections from hard bremsstrahlung, which in the present 
case corresponds to double hard bremsstrahlung. 

In the calculation of the contributions from the emission of two hard 
photons, it is convenient to decompose the phase space into same three re- 
gions discussed in 0: i) both hard photons hit the forward photon detec- 
tor, i.e., both are emitted within a narrow cone around the lepton beam 
($1,2 < ^o, $o *C 1); ii) one photon is tagged in the PD, while the other is 
collinear to the outgoing lepton (i9' 2 = Z(/c2,p') < $' ); and finally Hi) the 
second photon is emitted at large angles (i.e., outside the defined narrow 
cones) with respect to both incoming and outgoing lepton momenta. The 
last kinematic domain is denoted as the semi-collinear one. For the sake of 
simplicity, we shall always assume below that $g <C 1. 

The contribution from the kinematic region i) (both hard photons being 
tagged, and only the sum of their energies measured), has the form ||: 



i dv: 



77 



y dxdydz 



a 



L [ P%\z) + 2 1 -^- (\nz - ^ - 21b A 



+ 6(1 -z) + 
+ O (const) , 



4 



1-z 



1 - z In 2 z 



[l + z) 



In 



A 



(31) 



with 



P$\z) 



^21n(l-z)-lnz + ^ 



+ -(l+z) lnz-1 + z 



The contributions to the kinematic regions ii) and Hi) in the present 
case are slightly more complicated than for the purely leptonic measurement 
described in . Before discussing the contributions to the cross section from 
radiation almost collinear to the final state lepton, we shall therefore extend 
our treatment of the kinematics now and exhibit the necessary changes to 
the notation introduced in section 2. 
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To this end, let us recall the introduction (^) of the variable k. Again we 
define 

P-h 1 + cos^ 

K ■= irj = x *—2—- ^ 

Here x 2 is the energy of the second (i.e., non-collinear) photon that is not 
tagged in the PD, in units of the initial lepton energy, and d is its emission 
angle with respect to the incoming lepton. By similar reasoning as in the 
previous section, one can easily see that we presently have < k < z. 

On the other hand, the tagging of the collinear hard photon with energy 
in the PD corresponds to a reduction of the effective initial lepton energy 
by a factor of z (see eq. ^4|). The corresponding relations between measured 
quantities and invariant variables are simply obtained from eqs. (|T(J|)-(|r8|) 
via the simultaneous substitutions E e — > zE e , S — > zS, k — > k/z. Therefore, 
e.g., eqs. flTo] ) now read 

E' e = ( Z - K )(l-y) Ee + ^E p , 

Z — K 

(z-K) 2 (l-y)E e -xyzE p 

COSb' = -r- r . (66) 

(z - k) z (1 - y)E e + xyzE p 

In addition we define 

Q] = 2z P .p'=^ = X ^, 

Z — K Z — K 

i = -2zp ■ k 2 = 2x 2 zE 2 e {l - costf) , 
s = 2p ■ k 2 = 2x 2 E e E' e (1 - cos . 

The momentum transfer to the hadronic system and the true hadronic scaling 
variable are 

Ql = -(zp-p'-k 2 ) 2 = Q 2 -t-s, 

Xh = 91 = 91 (u) 

2P ■ (zp - p' - k 2 ) (z-k)yS' V ; 

The kinematic limit x 2 na,x ('&,'&') for the emission of the second hard photon 
is again obtained by requiring that the hadronic mass 

W 2 = (P + zp - p' - k 2 f = M 2 + ^—^Ql (35) 

be larger than the inelastic threshold, M 2 . More details can be found in 
appendix IB]. 
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The contribution to the radiative cross section from the semi-collinear 
region Hi) is obtained from the corresponding expression in the case of the 
e-method (see eqs. 3.7f of ref. M), 



1 d 3 cx 77 
_ in 

y dxdydz 



a p( t \ f dk2 a (Ql) rr, , U \ /qcn 
2 ( > o) / 77 n4 I 1 {zp } p } k 2 ), (36) 



with 

r 



GF 1 (z h ,Q*) + 



x ft [z 2 + 



>- K ) 2 (l-y) 2 ] S 2 -^-G 



+ 



(z- K )(l-£)(Q 2 -5)-z(4 2 -t) 



5 



F 2 (x h ,Ql) 



G = Qi-2St + Qt 



(37) 



and s,i,Qj,Q\, and as defined above. The angular part of the k 2 - 
integration is clearly restricted to the kinematic region Hi), i.e., the full solid 
angle with the exception of the separately treated cones around the incoming 
and outgoing lepton. 

Let us finally turn to the discussion of the kinematic region ii). As was 
already discussed in @, we expect that the contribution of this region to the 
observed cross section will depend on the experimental event selection, i.e., 
on the method of measurement of the scattered particles. We shall focus on 
the very same two possibilities. The first one is denoted as an exclusive (or 
bare) event selection, as only the scattered lepton is measured; the hard pho- 
ton that is emitted almost collinearly (i.e., within a small cone with opening 
angle 2$q around the momentum of the outgoing lepton) remains undetected 
or is not taken into account in the determination of the kinematic variables. 
The second, more realistic case (from an experimental point of view) is a 
calorimetric event selection, when only the sum of the energies of the out- 
going lepton and photon is actually measured if the photon momentum lies 
inside a small cone with opening angle 2d Q along the direction of the final 
lepton. 

First, in the case of an exclusive event selection, when only the scattered 
lepton is detected, we obtain for < 1 (similarly to eq. 3.3 of 



1 d 3 rr 77 

£ U °M,excl 

y dx dy dz 



cr 
4^ 



P{z,L ) 



^{x f ,y f ,Q 



C 2 



Co 



f) ' 



i + C 2 



L-l +(1-C) 



(3* 



(^We have of course kept the small but finite electron mass in the kinematic region of 
collinear radiation wherever necessary. See e.g. |15| for a very readable presentation. 
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where 



L 

no 

x f 



L' + 2 In 



no 



L' 



E 2 if 2 



77r 



+ 21nY(l), 



C 2 z(i-y) + [i-y(i-Q] 2 xye 

C [1 - y(l - 01 
[i-y(i-C)] 2 



o 



X . 



Vf 



(y 



i - m - o 



Q 2 



.f 



1 



C 2 

A 



l 



E p 



Y(l) ~ z(l-y)+xye 
and Co is the real solution in the interval (0, 1) of the cubic equation 

n 2 A T 



C 3 -£[l-y(l-Co)r 



[i-y(i-Co)]C 2 = o 



(39) 



(40) 



with A m = (M 2 — M 2 ) / S, see also appendix 0. We note in passing that the 
form of the leading logarithmic piece of (|38| ) agrees with the corresponding 



terms of the calculation performed in |L2 



Last, in the case of a calorimetric event selection, where only the sum of 
the energies of the outgoing lepton and photon is measured if the photon mo- 
mentum lies inside the small cone with opening angle 2d along the direction 
of the final lepton, the corresponding contribution reads 



1 d 3 



a 



77 



i,cal 



y dxdydz 



a 

4vr 2 



Cm ax 

p( Z ,L ) y dc 

o 



i + 
i-C 



(L' -l + 21nC) + (l-0 



(4-1) 2 In 



mi 

A 



2tt 2 



t. (41) 



The total contribution from QED radiative corrections is finally found by 
adding up (PSf), fl3"T|), (B§), and, depending on the chosen event selection, 



or fl41|). The reader may easily verify that the unphysical IR regularization 
parameter A cancels in the sum, as it should. 

It is important to note that the angle i9' plays only the role of an in- 
termediate regulator for the bare event selection and drops out in the final 
result. In the calorimetric case there are no large logarithmic contributions 
from final state radiation as long as $ does not become too small, since the 
mass singularity that is connected with the emission of the photon off the 
scattered lepton is canceled in accordance with the Kinoshita-Lee-Nauenberg 
theorem HTH1. For more details we refer the reader to PI. 
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4 Numerical Results 



In order to illustrate our results, we shall now present some numerical values 
obtained for the leading and next-to-leading radiative corrections. To facili- 
tate the comparison of the results for the S method with those for determi- 
nations of the kinematic variables based on a lepton-only measurement 0,0, 
we used as input 

E e = 27.5 GeV , E p = 820 GeV , $ = 0.5 mrad . (42) 

Unless stated otherwise, we chose the ALLM97 parameterization |T7| as 
structure function with R = 0, no cuts were applied to the photon phase 
space, and we assumed a calorimetric event selection. For the sake of sim- 
plicity we took a fixed representative angular resolution of i9q = 50 mrad 
for the electromagnetic calorimeter to separate nearby hits by an electron or 
positron and a hard photon, which is close to realistic for the HI detector at 
HERA. Also we disregard any effects due to the magnetic field bending the 
scattered charged lepton away from a collinear photon. 
Figure [I] compares the radiative correction 

dV 

5rc = tt 1 ( 43 ) 

with leading and next-to- leading logarithmic accuracy at x = 0.1 and x = 
10~ 4 and for a tagged energy of Epo = 5 GeV. The corresponding results 
for a tagged energy of Epp, = 20 GeV are shown in figure |j. The apparent 
cutoff at small y for small x is due to touching of the narrow cones defined 
by the solid angle covered by the photon detector and the cone around the 
final state lepton. 

At first sight the radiative corrections, being only of the order of a few per- 
cent, look rather small as compared to, e.g., a leptonic measurement. How- 
ever, this apparent suppression of the leading logarithmic part of the QED 
radiative corrections is easily traced back to the known weak dependence on 
initial state radiation of the determination of the kinematic variables (|2"B|)- 
(|25|) using the E method []6|,[T0f| . On the other hand, the pure next-to- leading 
logarithmic corrections are unsuppressed, as this "cancellation mechanism" 
does not work for non-collinear photon radiation. 

When we choose a bare event selection instead of a calorimetric one, 
the radiative corrections do become slightly larger. This is demonstrated in 
figure | that compares the radiative corrections for both selection schemes. 
The difference between the curves may be easily understood by noticing that 
in the calorimetric case the contribution from final state corrections to the 



13 



cross section is proportional to the relatively small logarithm ln$' , while it 
depends on the larger logarithm In E' e /m in the bare case. 

Next we shall study the influence of a photon energy cut on the radiative 
corrections. For simplicity, we assume an emission angle independent cut 
E™ x , which may be realized at large angles by rejecting events that show 
energy in the electromagnetic calorimeter sufficiently separated from the final 
state lepton, and for small angles by a cut on the variable 

S := Y> h + E' e (l + cos 9) - 2{E e - E PD ) , 

which is about twice the energy of a photon that is lost outside the PD. 
Figures f| and |^ illustrate the dependence of the radiative corrections on an 
energy cut on the semi-collinear (lost) photon. The influence of a rather loose 
cut of E 7 2 < 5 GeV is significant, especially at the larger value of x = 0.1, 
although the inclusive corrections were seen to be quite small. 

Finally, figure [| shows the dependence of the next-to-leading logarithmic 
corrections on the ratio R (see eq. p7|) , again for x — 0.1 and x = 10~ 4 and 
for a tagged energy of E PD = 5 GeV. As one would expect, one sees that 
only for large y there is a visible difference between the corrections calculated 
for (assumed constant) R = and R = 0.3 of the order of a permille. In 
the case of the purely leading logarithmic corrections the change would be 
much less than the order of the line width, so we omitted the corresponding 
lines. For this reason, a poorly known R(x, Q 2 ) as input to the calculation of 
the QED corrections will not have any significant effect on the extraction Q 
of R from the measured tagged photon cross section. Increasing the value 
of R up to, say, R = 1 for the smaller value of x would simply increase the 
difference with respect to the curve for R = but not lead to a qualitative 
change of our conclusion. 

5 Summary 

The £ method for the determination of the kinematic variables in deep in- 
elastic scattering has been reviewed in some detail from a theoretical point 
of view. We derived the relevant kinematics for the calculation of the hard 
photon emission contributions of the QED corrections to deep inelastic scat- 
tering for the HERA collider. As an application, we extended the formalism 
to radiative DIS events with a hard photon tagged in the forward photon 
detectors of the HI and ZEUS experiments. We have adapted the calcu- 
lations of the radiative corrections to these DIS events with a tagged 
photon for the £ method determination of kinematic variables. It turned 
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out that for a calorimetric measurement of the final state lepton the leading- 
logarithmic corrections are suppressed and thus quite small (of the order of 
5%), which is an intrinsic feature of the £ method. However, the typical 
size of the unsuppressed next-to-leading logarithmic contributions is of the 
same order of magnitude. However, for a bare event selection there are also 
significant contributions to the corrections already at the leading logarithmic 
level. 

The smallness of the leading QED corrections for the £ method suggests 
that the corrections are well under control, so that neglected higher order 
corrections will not play a significant role. It is also encouraging that the 
dependence of the corrections on the poorly known longitudinal structure 
function Fi (and thus R) is very small. 
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A Kinematic Limit of Hard Photon Emission 

This appendix is devoted to a discussion of the phase space limit for the 
emission of a hard photon in radiative DIS in the case of the S method. 
As we shall refer exclusively to the £ method, and to avoid cumbersome 
notation, we drop the index E from all kinematic variables. 

In the parameterization of the photon phase space following from @, the 
invariant mass of the hadronic system ( |IBD is 




+ 2 Xl E e E' e (l-cos$') 
M 2 + yS [1 - k - x] - 4£ e 2 



l2 1 — COS'i? 



1 + cos "d 



«[1-(1 -«)(!_„)] 




(44) 



where we heavily used the relations (|112|)-(|T5|) and 



1 + cos -d 



(45) 



Hi — X, 



7 



2 
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With the help of the following abbreviations 

M 2 -M 2 A 



a = 1 — x = 1 — x 



in 



X 



'E e 1 - costf 
E p 1 + cos $ 



v = \j \xy{\ — y) cosy? , (46) 
which satisfy the constraints 

< n < 1 , < A < oo , \v\ < 1, (47) 

the inequality P^ 2 > M 2 can be brought into a simpler form: 

- [y(l + A 2 ) + Xv] k - A 2 (l - y)K 2 > . (48) 

Obviously the physical range is < k < Ko, with 

y(l + \ 2 ) + \u ^[ y ^ + X 2 ) + Xv] 2 + Ay{l-y)^Y 2 
K0 ~~ 2X 2 (l-y) + W(l^y) • (49) 



Direct inspection of ([48]) shows that for y — > 1, Ko approaches the value 

k = K (y=l)= 1 - X -^ m <l. (50) 
1 + X z 

Note that the choice of $ and tp for the parameterization of the photon 
phase space is not well suited for analytic or semi-analytic calculations, as 
the treatment of the separation between the phase space regions ii) and in) 
in section 3 is rather cumbersome. 



B Parameterization of the Hard Photon Phase 
Space 

Instead of using the parameterization (f|) for the phase space of the hard 
photon, it is often convenient to trade the azimuthal angle p in favor of 
the angle *&' between the photon and the final state lepton. Introducing the 
variables n, r 2 , 

1 — cos $ 



1 — cos *&' 1 — cos 9 cos d — sin 9 sin d cos 09 

T2 ■= ^ = ^ -, (51) 
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the integration over the photon solid angle becomes 

Jdtl^ yd(cOS^)d^ = Jj{T 1 ,T 2 )e{V)dT 1 dT 2 , (52) 



with the Jacobian 

J(n,T 2 ) = 



8 



sin 6 sin d sin if \ y/j) 



where 

V = -±tt x t 2 - (r 2 + r 2 + r 2 2 - 2rr x - 2tt 2 - 2t x t 2 ) , (53) 
1 — cos 9 xye . . 

T := 2 = (l-K) 2 (l-y) + xye' ( } 

£7p 



e 



A factor of 2 has been taken into account in the Jacobian for the two-fold 
ambiguity in the azimuthal angle ip. Note that < k < 1 implies 

r>r<°>(s,y) = — . (55) 

1 — y + xye 

The range of integration for the variables T\ )2 for a given photon energy x^E e 
follows from the argument of the step function 0(T>). It is trivially obtained 
if r depends only on the kinematic variables x, y but not on the photon 
phase space variables, as is the case for a purely leptonic measurement of the 
kinematic variables 0, but it is more involved for the £ method and will 
therefore be discussed below. 

As the relevant variables for the parameterization of the phase space 



of the photon we choose r±, r 2 , and x 1 . Due to ([54]), r is related to x 



-7) 



aS K Xry (1 — Ti). Therefore, it appears to be reasonable to take T\ as the 

outermost integration and to determine the range of integration for the other 
two variables for each value of Ti, < T\ < 1. 

Obviously, the range where the argument T> ( |53|) of the G function above 
is positive is given by the interior of an ellipse as shown in figure [F[ The 
ellipse touches the r- and r 2 -axes at the value of Ti, and the lines r = 1 and 
t 2 — 1 at the value of (1 — n). 

The explicit form of the upper and lower boundaries of the ellipse are 

T±(n,T) = 7i(l - r) + r(l - n) ± 2y/ ri (l - r)r(l - n) . (56) 

On the other hand, remembering that for given x, y there is a lower limit 
(|55|) on r, we see that the part of this ellipse left to the dotted line is certainly 
unphysical. 
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Yet we have not made use of a lower limit on the hadronic mass that 
is due to the inelastic threshold or to a lower cut on the invariant hadronic 
mass. To this end, we express the hadronic mass (031) in terms of the r 12 : 



W 2 



M 2 + yS 1 - re - 



x 



1 - K 



M 2 +yS 

72 



1 - re 



x 



1 - re 



- 4x 7 EgTi + Ax^E e E' e T 2 

S T\ 



e 1 — r\ 



(57) 



+ 5 



l-n 



1 - K)(l-y) xy 



1 - K 



K . 



Requiring W 2 > M 2 leads to the inequality 



W — M 



where 



A 

B 

C 
D 



[Are 3 + Bk 2 + Ck + D] > 

(1 - re)(l - n) 

(l-2/)r 2 



;i-ri)j/ + 



n - 2(1 - y)r 2 



(5f 



.Tyr 2 - 2y(i - n) + (l - n)A m + 

(l-n) [(l-x)y-A m ] . 



[1 - y)r 2 - n 



Since A > 0, D > 0, the cubic equation Are 3 + i?re 2 + Ck + D = will have 
either one or three real solutions: always a negative one which is obviously 
unphysical, and two positive or complex conjugate ones. 

In case there are three real solutions, only the smaller one of the two 
positive solutions may be physical. This may be seen by studying the limit 
T\ — > 1, where re has to go to for any finite photon energy x 1 . 

For small values of r 2 , i.e., radiation collinear to the outgoing lepton, one 
always has three real solutions. This is easily verified by direct inspection 
of Q58D for A — ► 0: the negative solution goes to — oo essentially as —B/A, 
and the resulting quadratic equation in the limit r 2 = always has a positive 
discriminant. 

In the remaining cases when there is only one real (and negative) solution, 
the effective upper limit on re may be indirectly obtained from (|54[) and the 
maximum allowed value for function of Ti j2 , similarly to (|56|), 



re < re(ri,r 2 ) 



1 - T-H 



-(0) 



T 



(0) 



(59) 
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which is entirely of geometric origin. 

Let us now compare the above limits that were derived from the lower 
limit (or a cut) on the hadronic mass with the kinematic limits for collinear 
photon radiation, as given in [[H| . The case of collinear initial state radiation 
(ISR) is recovered by setting t\ — > 0, so that k — > x 7 , T 2 E' e — > xyE p /(l — x 7 ), 
and from (p7|) one immediately obtains 

x y < 1-x . (60) 

y 



This is consistent with the limit onz (= 1 — x 7 ) given in [12| up to terms of or- 
der O (A m ) , as these authors did not consider neither the inelastic threshold 
nor a cut on the hadronic mass. 

Collinear final state radiation (FSR) corresponds to taking the limit r 2 — ► 
0, Ti — ► t = xye/ {[1 — x 7 (l — t)] 2 (1 — y) + xye}, with r being implicitly 
defined. However, in order to be able to compare our limit with the one 
given in ref. ||12|| , we shall think of the collinear photon as taking away the 
fraction (1— £) of the outgoing lepton-photon system, while the lepton retains 
the fraction (. Hence, we have k — [(1 — ()/(}p', and we may write 

P ■ k = — - — P ■ p => = 1 — y , 1 — k 



C r i-«i-C 1-1/(1-0 ' 

(61) 



and thus 



. 1-C , xy 1-C 1 xt/[l-y(l-Q](l-C) 1 Q (ao , 
V ■ k = — - — p ■ p = — • -S = — — — • -S . 62 

Inserting these relations into (p!8[) we obtain 



W 2 



= M 2 + C \* [1 ® ]2 yS > M 2 ■ (63) 

FSR C 2 [l — 2/(1 — 0] 

This inequality leads to a lower limit Co, < Co < 1> that is the single real 
solution of a cubic equation, W 2 ((q) = M 2 . The kinematic limit that follows 



from (|63| ) is consistent with the limit for FSR given in [[12) , again up to terms 



of order O (A m J that have been neglected in their work. 

Finally, we should mention that the above considerations directly apply 
only to the case of conventional deep inelastic scattering with a single radiated 
photon. However, the case of DIS with a tagged collinear photon and a hard 
non-collinear second photon is easily recovered, provided we perform the 
simultaneous substitutions 

E e — ► %E e , S — ^ zS , Xy — ^ — , k — , s — , and A m 

z z z z 

(64) 

in the expressions given above. 
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Figure 1: Radiative corrections 5rc (EH) with leading and next-to- leading 
logarithmic accuracy at x = 0.1 and x = 10~ 4 and a tagged photon energy 
of 5 GeV. No cuts have been applied to the phase space of the second 
(semi-collinear) photon. 
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Figure 2: Radiative corrections <5rc ([EJ) with leading and next-to- leading 
logarithmic accuracy at x = 0.1 and x = 10~ 4 and a tagged photon energy of 
20 GeV. No cuts have been applied to the phase space of the semi-collinear 
photon. 
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Figure 3: Comparison of the radiative corrections for calorimetric ("calo") 
vs. bare measurement of the scattered lepton at x — 10~ 4 and at x — 0.1 for 
a tagged photon energy of 5 GeV. 
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Figure 4: Comparison of the cut dependence of the radiative corrections at 
x — 0.1 for a tagged photon energy of 5 GeV. A cut of < 5 GeV has 
been applied to the phase space of the semi-collinear photon. 
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Figure 5: Comparison of the cut dependence of the radiative corrections at 
x = 10~ 4 for a tagged photon energy of 5 GeV. A cut of < 5 GeV has 
been applied to the phase space of the semi-collinear photon. 
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Figure 6: Comparison of the R dependence of the NLL corrections at x = 0.1 
and x = 10~ 4 for a tagged photon energy of 5 GeV. No cut has been applied 
to the phase space of the semi-collinear photon. 



26 




Figure 7: Schematic view of the kinematic range for the integration over the 
angles r and r 2 for fixed t\. The allowed range from the change of variables 
( |5T| ) alone is given by the interior of the ellipse. In addition, for given x, y, 
the region left to the dotted line is unphysical, as r > r^°'(x,y) from (|55|) . 
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